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Abstract 

We study the closed universe recollapse conjecture for positively curved Friedmann- 
Robertson- Walker (FRW) models in the Jordan frame of the second order gravity 
theory. We analyse the late time evolution of the model with the methods of the dy- 
namical systems. We find that an initially expanding closed FRW universe, starting 
close to the Minkowski spacetime, may exhibit oscillatory behaviour. 

1 Introduction 

Since the discovery that general relativity (GR) can been derived from an action principle, 
many nonlinear Lagrangians have been used for the construction of a metric theory of 
gravity. The reasons for considering higher order generalisations of the Einstein-Hilbert 
action are multiple. Firstly, there is no a priori physical reason to restrict the gravitational 
Lagrangian to a linear function of the scalar curvature R. Secondly it is hoped that 
higher order Lagrangians would create a first approximation to an as yet unknown theory 
of quantum gravity. Thirdly one expects that, on approach to a spacetime singularity, 
curvature invariants of all orders ought to play an important dynamical role. Far from the 
singularity, when higher order corrections become negligible, one should recover general 
relativity. 

An other important motivation for considering generalised theories of gravity is the 
hope that these theories might exhibit better behaviour near singularities. General rel- 
ativity leads to singularities in the spacetimes of all known cosmological models with 
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ordinary matter. Higher order curvature corrections in the gravitational action may rec- 
tify the problem and lead to cosmological models free from such pathologies, at the cost of 
diverging from a FRW behavior at late times P . There is a resurgence of interest in such 
theories which naturally arise in string-theoretic considerations (cf. brane models with 
Gauss-Bonett terms |2j). These theories are considered as alternatives to general relativ- 
ity in an effort to explain the accelerating expansion of the universe (for a pedagogical 
review see 0). 

This talk is about the late evolution of a closed FRW universe in vacuum in the context 
of the quadratic gravity theory. Section 2 is a short exposition of the field equations 
obtained from a Lagrangian which is a function of the Ricci curvature. The conformal 
equivalence of these theories with the Einstein equations having a scalar field as a matter 
source is discussed. In section 3 we investigate the structure of the four- dimensional 
dynamical system describing the evolution of a closed FRW universe in quadratic gravity. 
We find that periodic solutions may prevent this universe from recollapse. 

2 Jordan versus Einstein frame 

We consider higher order gravity theories (HOG) in vacuum derived from Lagrangians of 
the form 

L = fiR)^, (1) 

where / is an arbitrary smooth function and R is the Ricci curvature. By varying L with 
respect to the metric tensor, we obtain the vacuum fourth-order field equations 

/' (R) R,, - i/ (R) g,, - V,VJ' (R) + g.^Df (R) = 0, (2) 

where □ = g^^V and a prime (') denotes differentiation with respect to R. It is well 
known (see ^), that under the conformal transformation 

Qt^u = f (R) 9t^u, (3) 

the field equations reduce to the Einstein field equations with a scalar field as a matter 
source, namely 

Gf,u = T^u (g, 0) , 

where 

T,, (g, 0) = d,<l)d,<P - ^g^, [{d<j)f - 2V (0)] , 

and 

0=y|ln/'(i?). (4) 



Assuming that (jH) can be solved for R to obtain a function R (0) , the potential of the 
scalar field is given by 

^ = ^W-/)- (5) 

The conformal equivalence between higher order gravity and general relativity is some- 
times expressed formally as "HOG = GR + scalar field". However this "equation" is an 
oversimplification of the picture. The two frames^ are mathematically equivalent, but 
physically they provide different theories. In the Jordan frame, gravity is described en- 
tirely by the metric g^y. In the Einstein frame, the scalar field exhibits a "non-metric" 
aspect of the gravitational interaction, reflecting the additional degree of freedom due to 
the higher order of the field equations in the Jordan frame. 

If matter fields, collectively denoted by are present in the Jordan frame, the corre- 
sponding field equations become 

fR^u - \fg,. - V.Vyf + g,yUr = T^y {g, v]/) . (6) 
The generalised Bianchi identities imply that 

{g,^) = 

Therefore, free particles move on geodesies in the Jordan frame. We pass to the Einstein 
frame by conformally transforming (jHl) to obtain 

G^y = T^y{g,^)+f^Ag,g,<il). (7) 

Now the Bianchi identities imply that 

Wf^y{g,g,<il)^0, WT^AaA)^^- 

Therefore, equations (|7j) are formally the Einstein equations, but this theory is not phys- 
ically equivalent to GR. Of course one could firstly conformally transform the vacuum 
Lagrangian into the Einstein frame to obtain 

L = (r - r''d^<pdy<i) - 2V (</))) , 

and then add the matter Lagrangian Lm (g, with minimal coupling. In that case, the 
stress-energy tensors of both fields (p and \& are separately conserved, 

= 0, V^T^. (g, 0) = 0. 

^The term frame denotes the set of dynamical variables used and is not associated to any coordi- 
nate reference frame. In the literature, the original set of variables is called the Jordan frame and the 
conformally transformed set is called the Einstein frame. 



This discussion raises the question: which is the physical metric? There is no universally 
acceptable answer to this question in the literature (see [HI EI for a thorough presentation 
of different views). 

There is another subtle issue, rarely discussed. The conformal transformation 
preserves the causal structure of the two frames only if /' (R) is non-negative on the 
whole spacetime. A weaker assumption is that /' has the same sign on a connected open 
subset of the spacetime. However, both assumptions rely on the knowledge of the function 
R, and this is possible only after the field equations have been solved. In any case, the 
conformal transformation ^ may fail to be regular at all points of the spacetime. 

3 Cosmological considerations 

Since it is easier to tackle a problem in the context of general relativity than using the 
fourth-order equations Q, the conformal equivalence theorem allows certain results valid 
in general relativity to be transferred in HOG. If a problem can been solved in the Jordan 
frame it should be interesting to compare this result with the solution of the same problem 
obtained in the Einstein frame. 

We investigate the closed-universe recollapse conjecture for the / {R) = R+[3R^ theory. 
This conjecture states roughly that a closed universe cannot expand for ever, provided 
that the matter content satisfies some energy conditions and has non-negative pressures. 
For homogeneous and isotropic spacetimes in the Einstein frame this conjecture was found 
to be true for vacuum models |H] and for models with a perfect fluid with equation of state 
p = (7 — l)p, 2/3 < 7 < 2 [S]. More precisely it was shown that, an initially expanding 
closed FEW universe, starting close to the Minkowski spacetime cannot avoid recollapse. 
For the R + (3R^ theory in the Einstein frame in vacuum, the state of the system is defined 
by ^a, d, 0, 0^ G and one has a problem of scalar-field cosmology with a potential (cf. 
equation (jS))) 




From the general field equations Q we obtain the trace equation 



R + 3HR+ —R = 0, 



and the — equation 



R 



R^ 



-RH - — - 
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in vacuum. Setting x = 1/a our dynamical system in vacuum is 



R = v, 
X = —xH, 

H = -R-2H'^ - kx^, 
6 



subject to the constraint 



;i + 2(3 R) (H^ + kx^) + 2(3 Hv = ^R\ (9) 

6 



Since we are interested only for the closed, k = +1, models, from now on we omit k from 
the formulas. 

The only equilibrium point is the origin (0, 0, 0, 0). Linearisation shows that the eigen- 
values of the Jacobian matrix at the origin are ±i/v^6^, 0, 0. For nonhyperbolic equilibria 
the Hartman-Grobman Theorem does not give any information regarding their stability 
and therefore, we shall try to find the normal form of the system. Let T be the matrix 
which transforms the linear part of the vector field into Jordan canonical form. We write 
(jS)) in vector notation as 

z = Az + F(z), (10) 

where A is the linear part of the vector field and F (0) = 0. Using the matrix T, we define 
new variables, {u, w, y, x) = x, by the equations 

yp'^' ^ = H = w + y, x = x, (11) 

or in vector notation z = Tx, so that (fTUj) becomes 

X = T-^ATx + T-^F (Tx) . 

Denoting the canonical form of A by -B we finally obtain the system 

x = Ex + f(x), (12) 

where f (x) := T~^F (Tx) . In components system ()12j] is 

1 

u = j=w, 

w = u — — 3wy, (13) 

v6p 



y = w — wy — 2y — x 
X = —wx — yx. 



Under the non-linear change of variables 
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y -^y — ^ V^«?/> (14) 

X ^ X + >/6(3ux, 

and keeping only terms up to second order, the system ()13|1 transforms to 

1 3 

u = ]=w yu, 

1 3 
w = , — u yw, 

y = \{u^ + w^)-2y^-x\ (15) 
X = —yx. 

Remark. The results of the normal form theory are valid only near the origin. Fur- 
thermore, the system (jH)) is not an arbitrary "free" four-dimensional system. In view of 
the — equation, the variables have to satisfy the constraint Q- The nonlinear coordi- 
nate transformation ()14|) is valid in a neighbourhood of the origin, small enough to ensure 
that the new variables still respect the constraint (jH)). 

Defining cylindrical coordinates {u = r cos 6,w = r sin 6,y = y,x = x) , we obtain 

3 

r = -ify^ 



6/3 

y = V-2y2_a;2, (16) 
X = —yx. 

We note that the 6 dependence of the vector field has been eliminated, so that we can 
study the system in the (r, y, x) space. The second equation of (fTB|) implies that the 
trajectory in the u — w plane spirals with angular velocity 1/a/6^. We write the first and 
fourth of (Uni) differential equation 

dr 3 r 
dx 2 a;' 

which has the general solution 



r = Ax^''^, A > 0. 



(17) 



We substitute (fTTj) into the third equation of (fT^ and we obtain the projection of the 
forth-dimensional system on the x — y plane, namely 



X = —yx, 

y = bx^ -2y^ -x^, b > 0. (18) 

Some general comments about the trajectories of the system (|T8|l follow by inspection. 
Firstly, (fT^ is invariant under the transformation t — > —t, y — > —y, which implies that 
all trajectories are symmetric with respect to the x axis. Secondly, the line a; = is 
invariant and therefore, a trajectory starting at the half plane x > remains there for 
all t > 0. Note that y is decreasing along the orbits in the strip < a; < 1/6 while x is 
decreasing in the first quadrant. On any orbit starting in the first quadrant with x < 1/6, 
y becomes zero at some time and the trajectory crosses vertically the x— axis. Once the 
trajectory enters the second quadrant, x increases and y decreases. System (fTHj) has two 
equilibrium points, the origin (0,0) and and (1/6,0) on the x— axis. Computation of 
the Jacobian matrix at the equilibrium (1/6,0) shows that this point is a center for the 
linearised system. Hopefully (fTK|l has a first integral, viz. 

A.( ^ 26 1 y2 

X X"^ X* 

In fact, writing (fTHjl as 

dy bx^ — 2y'^ — x^ 



dx —yx 

and setting y"^ = z, we obtain a linear differential equation for z which is easily integrable. 
The level curves of are the trajectories of the system. 

We shall show for the system (fTH|) that (i) there are no solution curves asymptotically 
approaching the origin (ii) there exist periodic solutions and (Hi) the basin of attraction 
of every periodic trajectory is the set y"^ + x^ — 2bx^ < 0. 

Proof. The function (p has a local isolated minimum at (1/6,0) and therefore its level 
curves near this point are closed. For (x, y) = C we have 

y^ = x^ (-l + 26x + Cx2) , 

which implies that —1 + 26x + Cx^ must be non-negative. It follows that for C > any 
orbit starting in the first quadrant satisfies 

x>^[~b + V^Tc'^ > 0, 

i.e., there are no solutions approaching (0,0) . For C G (—6^,0) an orbit of (fTHj) starting 
in the first quadrant crosses the at (—6 — y/W+C) /C and re-enters in the first 

quadrant crossing the at (—6 + V6^ + C) / C, i.e. it is a closed curve and represents 

a periodic solution. The curve corresponding to C = separates the phase space into 



two disjoint regions I and II. In region I every initially expanding universe eventually 
recoUapses. In region II, (C < 0), every trajectory corresponds to a periodic solution and 
we conclude that the basin of attraction of every periodic trajectory is the set + — 
2bx^ < 0. ■ 

The phase portrait is shown in Figure 1. 




Figure 1: Phase portrait of (fTHj) 



4 Discussion 

The periodic solutions of (fTH|) induce periodicity to the full four- dimensional system (fTB|l 
or (fT3j). In fact, if x (t) and y (t) are periodic solutions, then (fT7|) implies that the functions 

u(t)=r (t) cos , w = r{t) sin (^t) 

oscillate in the u — w plane with a periodic amplitude r{t). Obviously one cannot assign 
a physical meaning to the new variables (m, w, x, y) since the transformations (lllj) and 
(fT^ have "mixed" the original variables of (jH)) in a nontrivial way. However, the periodic 
character of the solutions of (fT3j) whatever the physical meaning of the variables be, has 
the following interpretation. Close to the equilibrium of the original system there 
exist periodic solutions for all variables. This implies that an initially expanding closed 
universe can avoid recollapse through an infinite sequence of successive expansions and 
contractions. This interesting result was not revealed in the Einstein frame jH] . Since the 
basin of attraction of all periodic trajectories of ()18p is an open subset of the phase space, 
there is enough room in the set of initial data of (jSj) which lead to an oscillating scale 



factor. It should be interesting to investigate the late time evolution of closed FRW models 
in the Jordan frame by including matter satisfying some plausible energy conditions. This 
issue will be presented elsewhere. 
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